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Anisotropic single-molecule magnets may be thought of as molecular switches, with possible ap- 
plications to molecular spintronics. In this paper, we consider current-induced switching in single- 
molecule junctions containing an anisotropic magnetic molecule. We assume that the carriers inter- 
act with the magnetic molecule through the exchange interaction and focus on the regime of high 
currents in which the molecular spin dynamics is slow compared to the time which the electrons 
spend on the molecule. In this limit, the molecular spin obeys a non-equilibrium Langevin equa- 
tion which takes the form of a generalized Landau-Lifshitz-Gilbert equation and which we derive 
microscopically by means of a non-equilibrium Born-Oppenheimer approximation. We exploit this 
Langevin equation to identify the relevant switching mechanisms and to derive the current-induced 
switching rates. As a byproduct, we also derive S-matrix expressions for the various torques entering 
into the Landau-Lifshitz-Gilbert equation which generalize previous expressions in the literature to 
non-equilibrium situations. 

PACS numbers: 73.63.-b,75.76.+j,85.75.-d 



I. INTRODUCTION 

In recent years, electronic transport through nanos- 
tructures has witnessed a shift toward molecular sys- 
tems. Several ingenious schemes for measuring transport 
through single molecules have been realized and experi- 
mental control over such systems is rapidly improving.'^' 
A prime difference between transport through single 
molecules or nanoelectromechanical systems (NEMS)^ as 
opposed to transport through more conventional nanos- 
tructures lies in the coupling of the electronic degrees 
of freedom responsible for transport to few well-defined 
collective modes of the molecule, with recent research 
focusing on effects of molecular vibrations (molecular 
nanoelectromechanics^pMl and local magnetic moments 
(molecular spintronics). ^SJIiSl 

The interesting property of transport setups based on 
single-molecule magnets is the possibility of combining 
some classical properties of macroscopic magnets with 
quantum features such as quantum tunneling. Work 
on molecular spintronics has focused on single molecule 
magnets such as Mni2 and transition metal complexes. 
Transport experiments with Mni2 concentrated on sig- 
natures of the magnetic excitations, as revealed by peaks 
in the differential conductancejISl and a spin-blockade 
mechanism. Research on transition metal complexes, 
based e.g. on Co, also addresses the Kondo effect. Re- 
lated phenomena have been discussed in molecular spin 
valves, which have been realized in setups with Cgo,'^ 
and more recently in TbPc2 setups coupled to nanotubes 
through supramolecular interactions 

In addition to their remarkable fundamental quan- 
tum transport properties, single molecule magnets are 
also appealing for their potential as memory cells in 
spintronics. In this context it is important to have re- 



liable mechanisms for writing and reading the stored in- 
formation. Specifically, it is essential to have protocols 
for manipulating and for detecting the orientation of the 
magnetic moment. To this end it is convenient to take 
advantage of the coupling between the spin of the elec- 
trons, which tunnel from the electrodes, and the localized 
magnetic moment of the molecule. 

Much of the existing literature, both on molecular na- 
noelectromechanics and molecular spintronics, assumes 
that the electrons reside on the molecule for times large 
compared to typical vibrational or magnetic precession 
periods. In this limit, it is often appropriate to treat 
the dynamics of the system within a rate- or Master 
equation in the exact eigenstate basis of the isolated 
molecule,'^'^^ describing also spin-transfer torques out-of 
equilibrium.'iSHII] ti^e context of nanoelectromechan- 
ics, there has recently been much interest in the oppo- 
site regime of adiabatic vibrational dynamics, in which 
the electronic processes are fast compared to the vibra- 
tional degrees of freedom, e.g., in the context of certain 
molecular s witc hes, NEMS near continuous mechan- 
ical instabilities,'2321] flexural modes of suspended carbon 
nanotubes and graphene,'^ or the cooling and amplifica- 
tion of mechanical motion by the backaction of conduc- 
tion electrons.'^ZHH 

The goal of the present work is to explore the latter 
limit in the context of magnetic molecules. We consider a 
generic model for the magnetic molecule, which includes 
an easy axis anisotropy sandwiched between two metal- 
lic (possibly polarized) electrodes at which a bias volt- 
age is applied.!^ We focus on the regime where the typ- 
ical time for dynamics of the molecular magnetic mo- 
ment is much larger than the dwell time of the elec- 
trons fiowing through the structure. Within this adia- 
batic regime it is possible to study the coupled electronic 
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transport and spin dynamics within a non-equilibrium 
Born-Oppenheimer approximation (NEBO) analogo us to 
the one adopted in NEMS in the equivalent regimeJ22122l 
Starting from a microscopic description, we can derive 
semiclassical equations of motion for the local magnetic 
moment that have the structure of generalized Landau- 
Lifshitz-Gilbert (LLG) equations. The latter have been 
the basis for several previous works in spintronics and 
nanomagnetism.l^SBSI 'We note that in previous works 
describing magnetic nanoparticles, LLG equations have 
been derived in a perturbative way assuming either that 
the coupling between the electronic spin and th e mag- 
netic moment of the nanoparticle is smalP^I^HH] g^^^ j qj. 
that tunneling between the leads and nanoparticle is 
weak.^'^ In contrast, our microscopic derivation re- 
lies entirely on the non-equilibrium Born Oppenheimer 
approximation which is valid in the high-current limit as 
described above. As a consequence our non-perturbative 
approach allows us to compute how the parameters of 
the LLG equation depend on the state of the molecular 
moment as well as on the applied bias and gate voltages. 
We mainly focus on two important features. First we an- 
alyze the magnetic molecule attached to (non-magnetic) 
metallic leads. In this case, switching of the molecular 
moment is induced by the fluctuating torque exerted by 
the current flow. In addition, we also investigate the 
renormalization of the switching barrier by the average 
torque caused by the charge carriers. Second, we con- 
sider that switching is dominated by a different mecha- 
nism for spin-polarized electrodes, namely by the spin- 
transfer torque exerted by the transport current. The 
latter is well known in the context of layered magnetic 
structures.l33E£l We also analyze the behavior of the elec- 
tronic current and we identify in this quantity the in- 
terplay between the spin fluctuations and the signatures 
of coherent transport, which are typical of the molecular 
devices. 

This paper is organized as follows. In Sec. |IT] we in- 
troduce our model of the single-molecule junction con- 
taining an anisotropic magnetic molecule. The Landau- 
Lifshitz-Gilbert equation describing the dynamics of the 
local moment of the molecule is derived within the NEBO 



approximation in Sec. Ill and related to scattering matrix 
theory in Sec. |IV[ This Langevin equation is explored in 
Sec. |V]S witching of the molecular moment is discussed in 
Sees. |VI| and |VII[ Section |VI| fo cuses on switching caused 
by fluctuations while Sec. |VIl| discusses situations when 
the switching is dominated by the spin-transfer torque. 
We conclude in Sec. IVIlH Some technical details are rel- 
egated to appendices. 



is coupled to left iV) and right (i?) leads at different 
chemical potentials. The spin § of the current-carrying 
electrons couples to a localized molecular spin M through 
exchange. Then, the full Hamiltonian 



a — ^ {Ha + Hda) 



a=L,R 



encompasses the Hamiltonians 



(1) 



(2) 



of the left {L) and right {R) leads, modeled as free- 
electron systems (creation operators c|, We will con- 
sider the possibility of spin-polarized leads, assuming a 
spin-dependent dispersion efc^^o-- The tunneling Hamilto- 



H. 



da 



E 



X(i„.ad^ + dUk^,<y) 



(3) 



describes the hybridization between the molecular orbital 
(with creation operator dl.) and the leads. The molecular 
Hamiltonian is given by 



geS-B-f Js-M + C/(M). (4) 



The potential experienced by the molecular spin in the 
absence of coupling to the external leads is J7(M) — gdM.- 
B — DM^. The uniaxial anisotropy of the molecule is 
parametrized through the anisotropy parameter D, with 
easy-axis anisotropy corresponding to _D > and easy- 
plane anisotropy to < 0. The coupling constant J 
denotes the strength of the exchange interaction between 
the molecular spin M and the electronic spins. 



(5) 



where (with j = x,y,z) are the Pauli matrices. For 
simplicity, we assume this exchange interaction to be 
isotropic. The energy of the molecular orbital eq can 
be tuned by a gate voltage and B represents a Zeeman 
fleld acting on the electronic and the localized spins with 
^-factors ge and gd, respectively. 



III. DESCRIPTION OF THE SPIN DYNAMICS 



II. MODEL 

We consider a minimal model of an anisotropic 
magnetic molecule embedded into a single-molecule 
junction.*^ We assume that transport through the 
molecule is dominated by a single molecular orbital which 



We now discuss this model in the limit of slow pre- 
cession of the magnetic moment, that is, many elec- 
trons are passing the molecule during a single preces- 
sional period of the molecular spin. In this limit, it is 
natural to approximate the molecular spin as a classi- 
cal variable whose dynamics can be described within a 
non-equilibrium Born-Oppenheimer approximation. The 
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resulting dynamics takes the form of a Langevin equa- 
tion of the Landau-Lifshitz-Gilbert type which we de- 
rive microscopically for our model. Specifically, the ex- 
change coupling between the current-carrying electrons 
and the molecular moment introduces additional torques 
and damping terms which enter into the Langevin equa- 
tion and which we will now discuss in detail. 



A. Semiclassical equation of motion of the 
molecular spin 

Our derivation starts from the Heisenberg equation of 
motion for the molecular spin, 



^£jik[Jsi +gdBi]Mk 



Lk 



+D J2 £zjk[M^Mk + MkM,], 



(6) 



where Sjik is the antisymmetric Levy-Civita tensor. 
Within the non-equilibrium Born-Oppenheimer approxi- 
mation, we can turn this into an equation of motion for 
the expectation value M(t) = (M(t)) of the localized 
spin. 



M = M X [-dMU{M) -Js + SB] , 



(7) 



with -dMU{M) = -gaB + 2DM,e^. Here, M = M(t) 
denotes the molecular spin averaged over a time interval 
large compared to the electronic time scales, but small 
compared to the precessional dynamics of the molecular 
spin itself. The corresponding time-averaged electronic 
spin s = s{t) can be expressed in terms of the electronic 
lesser Green's function 



G<At,t') = i{dl{t')d„{t)) 



of the molecular orbital as 



ih 



Sj7M(t) 



tr[G<{t,t)a^] 



(8) 



(9) 



It is important to note that due to the Born- 
Oppenheimer approximation, the lesser Green's function 
must be evaluated for a given time dependence of the 
molecular spin M(t). As a result, the average electronic 
spin s{t) depends on the molecular spin at earlier times. 
This will be considered in more detail in the next subsec- 
tion. The instantaneous contribution gives rise to a force 
acting on the molecular spin. Retardation effects pro- 
duce terms proportional to M, appearing in the equation 
of motion as Gilbert damping and a change in the gyro- 
magnetic ratio. Additionally, fluctuations of the electron 
spin give rise to a fluctuating Zeeman fleld 63 acting on 
the molecular spin. 



B. Electronic Green's function in the adiabatic 
limit 

We now turn to evaluate the electronic lesser Green's 
function, accounting for the slowly varying molecular 
spin M(t). We start by considering the corresponding 
retarded Green's function 



G«,,(t,t') = -^0(i-O(K(t),4'(O})- 



(10) 



Since the electrons are non-interacting, we can obtain 
from G^ at the end of the calculation. From now on, we 
set h — 1. The retarded Green's function satisfies the 
Dyson equatiorP^ 

-J2GtAt^t')rT,r^, - hit') = S{t-t')S, 
Here we introduce the self-energy 

a 

with 



iiu{t—t') 



(11) 



(12) 



s^,.(-) = E 



Wk, 



17] 



(13) 



accounting for the coupling to the (possibly spin- 
polarized) leads, see also Appendix |A] It is convenient 
to introduce an effective magnetic field experienced by 
the electrons given by 



b(0 = -(JM(t)+.geB) 



(14) 



Notice that even if we consider a constant external mag- 
netic field, the effective magnetic field is time dependent 
due to the explicit time dependence of the molecular spin 
M = M(t). 

In order to implement the Born-Oppenheimer approx- 
imation, it is convenient to rewrite the Dyson equation 
in the mixed (Wigner) representation defined by 

d{t,u;) = y"dTe'"^0(t + r/2,t-T/2), (15) 

for a general quantity 0{ti,t2) depending on two times 
with central and relative times defined as < = {ti +t2)/2 
and T — ti—t2. The non-equilibrium Born-Oppenheimer 
approximation can now be implemented by noting that 
the dependence on the central time t is slow. Thus, con- 
volutions such as C{ti,t2) = J dia Ci(ti, <3)C2(t3,i2) can 
be approximated in Wigner representation through 



exp 
C1G2 



j(* 



G1C2 



CidtC2 - dtCid^C2 



(16) 
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in next-to-leading order using the shorthand Ci = 

For our problem, to lowest order in the slow changes 
of M, we then obtain for the Dyson equation 



1 



G« (c. - £0 - 



crb) 



-^a^G^9tT.b-^atG^(l-9^S«) (17) 

where = G^{t,u!) denotes the Green's function in 
the Wigner representation. In the above equation and 
in what follows, the Green's functions, as well as the 
self-energy, are matrices in spin space with elements 
G^^,{t,uj) and I]^^,(cj), respectively. In the strictly adi- 
abatic limit we drop the terms proportional to derivatives 
with respect to the central time. To this order we obtain 



G^(M, cj) = [cj - eo - S^(w) -a-h] 



(18) 



In next-to-leading order in the Born-Oppenheimer ap- 
proximation, we keep the time derivatives with respect 
to central time to linear order. Equation (18 1 implies 
[dtGQ']^^ = — cr • b. Accordingly, by differentiating 
[cli — eo — S]^(a;) — cr • b] = 1 with respect to time 
and multiplying the resulting equation with Gq one ob- 
tains the identity dtGQ — G^cr ■ hG^. Then the Dyson 
equation yields 



G' 



dujGnfT ■ i>Gn — Gn C 



hd^G^ 



(19) 



The lesser Green's function can now be deduced from 
the relation G< = G^' (E)J:< (g) G^,'3s! where (g) denotes in- 
tegration over internal time arguments and G^ — [G^]^. 
The lesser self-energy depends only on time differences. 



(20) 



Here, we introduced Ta.ai'^) ~ ^2Im[I]^^(a;)] as well 
as the Fermi functions /q(w) = 1/(1 -I- e'^°("~'^°)) with 
a = L,R. We obtain after straightforward algebra 



G" 



Gq 



[d^G<o- ■ hG^ - G<a- ■ hd^G^ 
d^G^a- ■ bG< - G^a- ■ hd^G<) . (21) 



Here we used Gq = Gq^Y^^Gq and suppressed the 
arguments of the frozen Green's functions, Gq' — 
Go«^^'<(M,c.). 



The first term in ( 22 ) contains the average electron spin 

duj 



s°(M) 



2n 



tr[GS^f7'i 



(23) 



in the strictly adiabatic limit. The correction due to re- 
tardation effects associated with the slow dynamics of the 
molecular spin are captured by the matrix 7(M), 

ynM)=^ / g{tr[a.G>a'=G<a'] 

+ tr[9^G^ {<j''G<a' - a'G<<7>')]}, (24) 

where we have integrated by parts and used the greater 
Green's function 



G>„,{t,t') = -^{d.{t)dl{t')), 



(25) 



with the relation — = G^ — G^ . It is appropriate 
to spHt this matrix into 7^^' =^ 7.^^ + 7^^' with the shorthand 



Is.a 



(7"=±7'=')/2, see Eqs. ( |A10| and ( |A11| . As we will 
see, the symmetric part 7^*^ describes Gilbert damping 
of the molecular spin, induced by the coupling to the 
electrons while the antisymmetric part 7^^^ will induce a 
coupling renormalization. 

Due to the stochastic nature of the current flow 
through the magnetic molecule (as reflected in thermal as 
well as shot noise of the current), the electronic spin will 
also fluctuate, giving rise to a fluctuating torque 53{t) 
acting on the molecular spin. Using Wick's theorem we 
obtain for the symmetrized correlator 

{5Bk{i)5B^(t')) = ^tr[a'=G>(t,i')^'"G<(t',t)], (26) 

of the electron spin. In the Born-Oppenheimer approx- 
imation, the fluctuations of the spin, as given by Eq. 



(26 1, can be evaluated using the Green's function Gg ' 
to lowest order in b. Thus, the fluctuating Zeeman 
field 5B has the symmetrized correlator {5Bk{t)5Bi{t')) — 
b^\M)5{t^t') with 



D^\M.) 



J2 

4 



2tx 



tr[a^G>a'G<] 



(27) 



Note that in the Born-Oppenheimer limit, we can neglect 
any frequency dependence of this correlation function on 
the time scales of the molecular spin, so that the fiuctu- 
ating Zeeman field can be taken as locally correlated in 
time. 



C. Electron spin 



D. Landau-Lifshitz-Gilbert equation 



We can now employ this result for the electronic 
Green's function and evaluate the electron spin. Sub- 
stituting Eq. ([21]) into Eq. ([9|, we find 



Js(M) ~ Js"(M) 7(M)M. 



(22) 



Substituting the expression for the electronic spin ( 22 1 
into the equation of motion ([7| we obtain a Langevin 
equation of the Landau-Lifshitz-Gilbert type, 

M = Mx [-SmC/- Js°-7.M-7aM + 5Bl . (28) 
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Note that, unlike in simple versions of a Landau-Lifshitz- 
Gilbert equation, the effective exchange field as well 
as the coefficient matrices 7^ and 7^ still depend on the 
molecular spin M itself. We can simplify this equation 
by introducing the vector 



CkiM) 



(29) 



Using that the length of M is conserved, it follows that 
the antisymmetric part of 7 merely renormalizes the pre- 
cession frequency by an overall prefactor 



a(M) 



1 



1 + C-M 



(30) 



This yields the simplified Landau-Lifshitz-Gilbert equa- 
tion 



M = aM X 



(31) 



which we will analyze further in the subsequent sections. 

When coupled to spin-polarized leads and when 
a finite bias voltage is applied, the torque can be 
non-conservative, yielding the so-called spin-transfer 
torque.*^ Also the eigenvalues of 7^ can become overall 
negative, providing another mechanism of energy transfer 
from the electrons to the localized spin. 

It is interesting to compare these results with those 
for the related problem of charge carriers interacting 
with a slow vibrational degree of freedom in a NEMS. 
In both cases, the dynamics of the slow collective de- 
gree of freedom can be described in terms of a Langevin 
equation .f22l2l] Since the stochastic spin dynamics is ef- 
fectively two-dimensional, it generically exhibits simi- 
lar phenomena as NEMS with more than one vibra- 
tional mode.l^^^ Specifically, this includes the non- 
conservative nature of the average force in general non- 
equilibrium situations as well as the presence of the anti- 
symmetric contribution to the velocity-dependent force. 
The latter Berry phase contribution'*^ acts, however, in 
different ways in the two cases, owing to the different 
orders of the Langevin equation. In the vibrational con- 
text, this term gives rise to an effective Lorentz force, 
while it merely renormalizes the precession frequency in 
the context of the magnetic molecule. 



IV. RELATION TO SCATTERING MATRIX 
THEORY 

Before proceeding with analyzing the Landau-Lifshitz- 
Gilbert equation (31) in more detail, we pause to pro- 



vide S-matrix expressions for the various entries into this 
equation. It has already been noted in a series of works 
by Brataas et aZ.I^SlESl that the coefficients in the Landau- 
Lifshitz-Gilbert equation in lead-ferromagnet-lead struc- 
tures can be expressed in terms of the scattering matrix of 
the structure, resulting in expressions for Gilbert damp- 
ing and the fluctuating torque in thermal equilibrium and 



for current-induced spin-transfer torques within linear re- 
sponse theory. Here we will provide S-matrix expressions 
which remain valid in general out-of-equilibrium situa- 
tions and which include the exchange field and the preces- 
sion renormalization in addition to the Gilbert damping 
with the only assumption that the precessional frequency 
of the localized magnetic moment is slow compared to the 
electronic time scales. Our discussion here closely follows 
recent work on current-induced forces in nanoelectrome- 
chanical systems.'^^'^ 

For adiabatic parameter variations, the full 5-matrix 
of mesoscopic conductors can be expressed in the Wigner 
representation as S{t,u}) = 1 - 2m{WG^W'^)it,uj). Ex- 
panding S to linear order in the velociti es M of the adi- 
abatic variables, Moskalets and Buttikeif^SlsiJ introduced 
an ^-matrix through 



S{t,ij) ~ S{M{t),uj) + M{t) ■ A{M{t),uj). 



(32) 



For the model considered here, the frozen 5-matrix is 
readily related to the frozen retarded Green's function 
G^(M,w) througlJ33 



S'(M,w) = 1 - 2TriWG^(M,uj)W^ 
while the ^-matrix is given by 
J 



(33) 



TT — 

2 



(34) 



The average electronic spin s^(M) can be written in 
terms of the frozen S-matrix ( 33 1 by expressing the lesser 



Green's function Gq 
energy Yi^{ui) 



G^^G^ in terms of the self- 



JTT^Q, faW^HaW with IIq, a projector 
on lead a. Using the identity 2mW'^W = (G^)"* - 



(G" 



we then find 



jsim = 



2m \ dMk 



(35) 



for the average electronic spin. Here the trace "Tr" acts 
in lead-channel space. 



The S-matrix expression ( 35 1 allows us to make some 



general statements about the average torque acting on 
the molecular spin. In particular, we can evaluate the 
curl of the average torque. 



5(J,s°) diJsf) 



dM, 



dMk 



de 



/aTr H 



dS'< dS 



(36) 



In thermal equilibrium, Eq. (36 1 can be turned into a 



trace over a commutator of finite-dimensional matrices 
due to the relations — f, H^ = 1, and unitar- 
ity S'tS' = 1. This implies that disc's) /dMk = so 
that there is no spin-transfer torque. In general out-of- 
equilibrium situations, the curl will be nonzero, giving 
rise to finite spin-transfer torque. 
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Figure 1: (Color online) Illustration of the setup. The elec- 
tronic level splits due to the effective magnetic field, given by 
b = b{&). The number of levels in the current window de- 
pends on b, the applied bias voltage eV — fiL — 1-i.r and the 
gate voltage eVg = (/ii, -|- aih)/2 — e. It changes from zero 
{eV/2 > eVg ± 6) to one {eVg + b > eV/2 > eVg - b) to two 
{eVf, ±b> eV/2) [assuming eV, eV^ > 0]. 



Similar to the average spin, we can also express the 



variance of the fluctuating Zeeman field ( 27 1 in terms of 
the frozen S'-matrix, 



7)'='(M)=5:|^/„(i-/„o 



dS 



dS 
dM, 



(37) 



By going to a basis in which D is diagonal and using 
Hq = n^, we find that I? is a positive definite matrix. 

To express the velocity-dependent forces in terms of 
the scattering matrix in general non-equilibrium situa- 
tions, we need to go beyond the frozen scattering ma- 
trix S and include the A matrix introduced above. The 
Gilbert-damping coefficients appearing in the Langevin 
equation (|31|) can then be written as 



,kl 



duj 



/aTr 



\dMi 



A, -A 



dS^ dS 
dMi 

dS 



(38) 



The eigenvalues of the first line are strictly positive while 
the sign of the second line is not fixed, giving rise to the 
possibility of overall negative Gilbert damping. Note that 
the second line is a pure non-equilibrium contribution. 
This ca n be seen by using unitarity of S as well as S, 
iniplyin^nmiAS't-HS'At = {i/2)[dMSd^S'< -d^SdmS^. 
With this preparation, it is now easy to ascertain that in 
equilibrium damping and fiuctuations are related by the 
fiuctuation-dissipation theorem, Z)'^' = 2Tj^''. 

Similarly, we express the antisymmetric part of 7^^' as 



dA,_dAi 
dMi dMi 




which causes a renormalization of the precession fre- 
quency, as discussed above. 
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Figure 2: (Color online) Component of the average current- 
induced torque in a uniaxial situation for unpolarized leads. 
Jsz{M.) is plotted as a function of the applied bias voltage 
for different gate voltages eVg. As discussed in the text, Jsz 
changes when the number of levels in the current window 
varies at Vg ± eV/2 = ±6 (see also Fig. [T|. The plots are 
obtained at zero temperature at the potential minimum i9 = 
for the values JM/2 = 0.2, F = 0.1, geB/2 = 0.002 and 
<7e = gd- The precession frequency in the absence of coupling 
to electron spin and magnetic field is = 2DM = 0.01. All 
energies are measured in units of the barrier height without 
magnetic field DM^ . 



V. MOLECULAR SWITCHES WITH AXIAL 
SYMMETRY 

From now on we specify to the case of axial symmetry, 
where both the magnetic field and the polarization of the 
leads point along the anisotropy axis. In this section, we 
will derive explicit expressions for the current-induced 
forces, including their dependence on the molecular spin 
M. 

We first consider the average torque which is deter- 
mined by the average electronic spin. Given that there 
are two basic vectors in the problem, namely e^, and 
M = M/M, the spin can be decomposed as 

s°(M) = sm(M)M + s,(M)e, + st(M)(e, x M). (40) 

Hence, the average torque exerted on the molecular spin 
by the conduction electrons is 



-am X Js" 



-aM X 



: + St (e^ X M) (41) 



which is obtained by inserting Eq. ( 40 1 into the Landau- 
Lifshitz-Gilbert equation (31 1. The first term inside the 



bracket can be derived from a potential, since its curl 
vanishes. This becomes more evident from the explicit 
expressions below using that the M-dependence of the 
coefficients stems from the effective magnetic field b ex- 
perienced by the electrons and that the length of M is 
conserved. This contribution modifies the precession fre- 
quency around the z-axis. In contrast, the second term 



on the right hand side of Eq. ( |41| has a non-vanishing 
curl. 



Vm X 



St (e^ X M) 



7^ 0, so that St introduces a 

non- conservative torque, providing the possibility of en- 
ergy exchange between the conduction electrons and the 
molecule. 
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Figure 3: (Color online) Damping and fluctuations as a function of the applied bias voltage, (a) The damping coefficient 7s, i 
is plotted for three different gate voltages. It is peaked when the number of levels in the current region changes from zero to 
one to two at ± eV/2 = ±6 (see also Fig. (b) At these points the fluctuation kernel Di increase steplike. The level 
broadening results from the interaction with the leads encapsulated in F. (c) The effective temperature T^s — D-^/ {2'ys,i) is 
shown as a function of the bias voltage. The plots are obtained at the potential minimum •& = Q with the same parameters as 
in Fig. [2] All energies are measured in units of the barrier height without magnetic fleld DM^. 



Concrete expressions for these contributions to the 
current-induced torque can be obtained from 



^ I g[G<(M,a.) 



s.{M.)=- 

+ G<(M,c.)(^+Re[I]f] 



(42) 
(43) 



as derived by substituting Eq. ( |A6| into Eq. (23 1 and 
taking into account possibly spin-polarized leads with the 
notation Y.^ ^^^^-^ = \Y^^ ± S^^]/2 for the self-energies. 

These general expressions simplify significantly for un- 
polarized leads, which corresponds to ^ = 0. Indeed, 
one finds that then Gf and Gf vanish, see Eq. ( A7 1. This 
implies in particular that the component st of the aver- 
age torque vanishes. The remaining conservative contri- 



bution Sz is then found to be s. 



arctan 





7r\ 


r/2 


+ 2) 



with 



9eB 

2b ■ 



(44) 



Here we assume the limit of zero temperature and intro- 
duce the shorthand e = eo + Re[Ec]- 

It is instructive to study the dependence of the average 
torque on bias and gate voltage. Notice that, due to the 
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Figure 4: (Color online) Potential experienced by the molec- 
ular spin in the absence of coupling to the electrons. f/(i?) is 
shown for different magnetic flelds. All energies are measured 
in units of the barrier height without magnetic held DM^ . 



effective magnetic field acting on the electron spin, the 
electronic level splits, see Fig. [T] The average torque is 
finite when just one level, corresponding to e.g. spin- 
up electrons, is occupied. In contrast, for sufficiently 
high bias voltages both spin-up and spin-down electrons 
participate in the transport so that no net electron spin 
acts on the molecule. This is illustrated in Fig. |2] where 
the average electronic spin on the molecule is plotted as a 
function of the applied bias voltage eV for three different 
values of the molecular level Cq (as tunable by the gate 
voltage eVg). 

For Gilbert damping and the fluctuating torque, we re- 
strict ourselves to unpolarized leads. This choice is mo- 
tivated by the fact that switching of the molecular spin 
(as discussed in the next two sections) is dominated by 
the average torque for polarized leads (and thus weakly 
affected by higher orders in the adiabatic expansion) and 
by the fluctuating force for unpolarized leads. (We men- 
tion in passing that expressions for Gilbert damping and 
fluctuating force for polarized leads can be readily de- 
rived but are rather cumbersome.) 

For unpolarized leads, we can split the Gilbert damp- 
ing tensor into one part proportional to the unit matrix 
and another proportional to a projector onto the z-axis, 

M X = 7s, iM X M + 7s,2(M • ez)M x e^, (45) 

where 7^ 1 and 7^ 2 are scalars. The first term in Eq. 
(45 1 tends to (anti-)align the molecular spin with the 



anisotropy axis while the second modifies the precession 
frequency. 

The coefficients 7^1 and 7^ 2 are calculated by inserting 
Go and Gq from Eq. ([AT]) into Eq. (|A12|, resulting in 



7a.l(M) 

and 

7s,2(M) = 



duj -k ^ 



e±&)2 + (r/2)2 



(46) 



du 



E 



{9eBJfTj:p{-d^fp){uj-ef 



n± [(a.-e±5)2 + (r/2)2] 



(47) 
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Figure 5: (Color online) Switching dynamics for unpolarized leads, (a) V is plotted for different bias and gate voltages, (b) 
The switching rate l/ri2 is plotted as a function of the applied bias voltage for different gate voltages, (c) The ratio between 
the switching rates l/r2i and l/ri2 is shown. The plots are obtained with the same parameters as in Fig. [5] All energies are 
measured in units of the barrier height without magnetic field DM^. 



The damping coefficient is peaked when the number of 
levels between fi^ and fiR changes and thus vanishes at 
large voltages when both levels are in the transport win- 
dow. We illustrate this dependence of 7^1 on gate and 
bias voltage in Fig.|3] The prefactor a in Eq. (31 1 is cal- 
culated in the same way as the damping coefficients, and 
the resulting expression is relegated to the appendix, see 
Eq. (|AT3|. 



for the variance of the fluctuating Zeeman field. Eq. ( 27 1 
which becomes Z)*'''(M) 



We close this section with the corresponding expression 
;uating Zeeman field. E 
Di{M)Ski + D2{M)bkbi, where 



^2(M) = J2 



27rf-n±[(^-£±&)' + (r/2)2] 



du! ^-v 

du x ^ 



(48) 



for unpolarized leads. As illustrated in Fig. [3j the 
strength of the fluctuations changes with the number of 
electronic levels in the transport window and saturates 
at high bias voltages when both levels lie within. 



VI. FLUCTUATION-INDUCED SWITCHING 
OF THE MOLECULAR MOMENT FOR 
UNPOLARIZED LEADS. 

We now apply our results to discuss the switching dy- 
namics for unpolarized leads, fn the absence of coupling 
to the electrons the molecular spin moves in the poten- 
tial U = gdBMz — DM^. For sufficiently small magnetic 
flelds, two minima are present, corresponding to parallel 
and antiparallel alignment of the spin to the magnetic 
fleld, see Fig. |4] 

Assume that the molecular spin is initially aligned par- 
allel to the magnetic fleld. Due to the interaction with the 
electrons the molecular spin fluctuates about this initial 
state, causing spin flips at a certain rate which we cal- 
culate in this section. Clearly, these fluctuations depend 
on temperature and applied bias voltage. If the system 
is in thermal equilibrium, this is a standard problem.'^ 



Our approach allows us to extend these standard results 
to out-of-equilibrium situations in the presence of a bias 
voltage in addition to flnite temperature. We also demon- 
strate that the orientation of the molecular spin can be 
read out by tracking the current through the molecule. 



A. Fokker-Planck equation 

Our approach is based on an equivalent Fokker-Planck 
formulation of the Langevin dynamics of the molecular 
spin. We flrst rewrite the Langevin equation (31 1 for un- 



polarized leads. Describing the orientation of the molec- 
ular spin in terms of a polar angle t?, measured relative 
to the applied magnetic fleld, and an azimuthal angle ip, 



and noting that M/A/ 
Langevin equation 



•de§ + ipsiwde^, we flnd the 



i} — a [M^s,i sin 1? — 5B^] 
sin-d tf — a \—d{)U/M — Js^ sini? + SB^ 
-A/(7,,i-H7,,2sin2 ^) 



(49) 



where the noise correlator is given in polar coordi- 
nates by i)'^'^ = D^'P = Di and D^^ ^ Di + 
[{geB /2) cos^ ain^f D2, with A deflned in Eq. (|48|. 

Following standard procedures,!^ we now derive the 
corresponding Fokker-Planck equation for the probability 
distribution P(M, t) of the magnetization vector M at 
time t. In the uniaxial situation under consideration, this 
probability distribution is independent of ip and depends 
on the angle only. As outlined in Appendix [B] for the 
convenience of the reader, we then obtain the Fokker- 
Planck equation 

dtP{i},t) = ^c»^sin?9e-^('')a^e^W/3(i9)P(?9,t). 
smw 

(50) 

This equation has the stationary solution P(i?)stat 
exp[— y(i?)]//3. Here we have introduced 



(51) 
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Figure 6: (Color online) Current-voltage characteristics, (a) The differential conductance dl/dV is peaked when the number 
of levels in the current window changes at Vg ± eV/2 = ±&. (b) Obviously, the current changes when the number of level with 
energy between the chemical potentials of the two leads changes, (c) The level splitting 26(i9) depends on the orientation of 
the molecular spin. The relative difference of the current evaluated at the two potential minima is plotted as a function of the 
bias voltage, [(a) and (b) are evaluated at i9 = 0. The plots are obtained with the same parameters as in Fig. |2] and all energies 
are measured in units of the barrier height without magnetic field DM^.] 



and 

VW= rM''-^.^'''^'^'^^''. (52) 
J Di/{2js,i) 

As long as the anisotropy is sufficiently large, U (??) has 
a minimum Ui at — — 0, another minimum U2 at 
"i? = '(?2 = TT and a maximum Um at < < tt. We 
assume that this holds also for V{-d) and visualize the 
dependence of ¥{1}) on gate and bias voltage in Fig. [s] 
One clearly sees that the difference between the values 
of V at the minima and the maximum decreases with 
increasing bias voltage, as one expects from the behavior 
of fluctuations and damping, cf. Fig. [3] 

Note that in equilibrium the ratio £'1/(275.1) = T, 
as dictated by the fluctuation-dissipation theorem. For 
zero temperature but finite bias voltages V it is some- 
times instructive to interpret this ratio as an effective 
temperature in each potential well, Tcff ~ Di/{2^s,i) (as 
done for instance in Refs. 22 23 36), see Fig.|3] Generally 
however, both coefficients, Di and 7^1 are angle depen- 
dent and non trivial functions of voltage, as we have seen 
explicitly above. 

We calculate how long the molecular spin remains on 
one half of the Bloch sphere. The mean time between 
passing the energy barrier r^j due to the interaction with 
the electrons is then found by a standard procedure. 



We consider an adjoint equation to Eq. (50 1 



m 

sini9 



e^(^)5^e-^(^) sinm n^i^) = -1/2, (53) 



with an absorbing boundary condition Tij{-d„i) = in or- 
der to get the mean first passage time, as briefiy outlined 
in Appendix [B] The factor 1/2 takes into account that it 
is equally likely to go to z? ^ at = i^m- Solving the 
equation yields 



= 2 



(54) 



for passing from 1? < 'd^n to z? > i?„i and an analogous 
expression for the opposite process. 



When the potential minima are well separated and the 
fluctuations are small, we can give an analytical expres- 
sion for the switching rate. In this limit, the integrals 
in ( |54[ ) can be evaluated by saddle-point integration (cp. 
Ret. 1521 for the situation in which the coefficients do not 
depend on d), yielding 



1 



1 



27r 



sin^™^|f"(^™)|^(T?,)F"(^,)e-[''^'''"^-''^'''^l. 

(55) 

Hence, the rate depends exponentially on the difference 
between V evaluated at its maximum and minimum, re- 
spectively, so that it can be tuned by varying bias volt- 
age and gate potential. The general behavior of 1 /r^- , as 



given by Eq. (54), is shown in Fig. [5] for typical values as 
a function of gate and bias voltages. We discussed above 
that the ffuctuations increase with the number of levels 
in the current window. This is also reffected in the ffuc- 
tuation induced transition rates which increase with the 
bias voltage accordingly. 

B. Current 

The current through lead a is given by the change of 
the number of particles in the lead times the electronic 



charge, = -ie{ H,J2k. 
batic limit this becomed^Sl 



2tt 



In the adia- 
(56) 

where = Sq^c for unpolarized leads. Noting that 
II = ~Ir — 1/2 and assuming symmetric coupling to 
the leads, Tl = = T/2, we obtain, by i nsert ing t he ex - 
pressions for the Green's functions [Eqs. (A4| and (A61] 
and the self-energies [Eqs. ( |Al| and ( pO| ] after straight- 
forward algebra 



e 

An 



arctan 



I r/2 



-{L^ R) 



(57) 
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Figure 7: (Color online) Sketch of the effect of polarized leads 
inducing spin-torque-transfer (indicated by the magenta ar- 
row) on the molecular moment. Depending on the polariza- 
tion and the current, the spin-torque tends to align the molec- 
ular moment either (a) along or (b) against the magnetic field. 



which is valid at zero temperature. As discussed above, 
the electronic level splits due to the interaction with the 
effective magnetic field b, defined in Eq. p4| . When this 
level splitting is larger than the level broadening F, the 
current increases as the number of levels in the trans- 
port window increases, see Fig. [T] This is reflected in 
peaks of the differential conductance dl/dV as a func- 
tion of gate and bias voltage. Note that the splitting of 
the electronic levels and thus the number of levels in the 
transport window depends on the molecular spin orien- 
tation since b = b{{}). As a consequence, the current is 
also a function of d. In principle, this allows one to read 
out the molecular switch via current measurements, see 
Fig.|6] 



VII. SPIN-TORQUE-INDUCED SWITCHING 
WITH POLARIZED LEADS 



The switching mechanism discussed in the previous 
section originates in fluctuations of the molecular mag- 
netic moment, introduced by the coupling to the itinerant 
electrons. In Sec. |IV| we have seen that the presence of 
polarized leads opens the possibility of negative Gilbert 
damping which could favor the switching of the molec- 
ular spin. This mechanism strongly depends on the de- 
tails of the system, like the value of the mean chemical 
potential fj, and the applied bias voltage. However, for 
spin-polarized leads, switching of the molecular moment 
under general non-equilibrium conditions will typically 
be dominated by a different mechanism which is driven 
by the non-conservative (or spin-transfer) torque exerted 
by the coupling to the current carrying electrons. The 
generic effect of the spin-torque in the dynamics of M 
has been reviewed in Ref. HSl This term appears already 
in leading order of the Born-Oppenheimer approxima- 
tion in which Gilbert damping and fluctuations can be 
neglected. 



Figure 8: (Color online) Motion of the molecular moment in 
the presence of spin-polarized leads. For negative bias voltage 
eV — —0.4, the magnetic moment is driven toward the posi- 
tive z axis (red curve) while inverting the voltage eV — > —eV 
causes a flip of the magnetic moment (blue curve). (We con- 
sider = F;^ = 0.1, F]: = r], = and eVg = 0; the other 
parameters are the same as in Fig. [2]) 



In this section we focus on this spin-torque —Jstle^ x 
M), see Eq. (41), in the Landau-Lifshitz-Gilbert equa- 
tion (31 1, where St is given by Eq. (43). We analyze un- 



der which microscopic conditions it is expected to drive 
switching in our molecular setup. In the present case it 
is clear that it moves the vector M along the azimuthal 
direction, tending to align it along the magnetic fleld. 
Thus, given a tilted molecular magnetic moment M pre- 
cessing around the magnetic field, for St < the spin 
torque induces a spiral trajectory moving M toward or- 
bits of smaller radius around the magnetic field. Instead, 
for Sf > it induces orbits of larger radius enabling the 
switching to the opposite hemisphere, with M tending to 
align opposite to the external magnetic field. 

In our model, the behavior of st can be rather easily 
analyzed in the limit of completely polarized leads, e.g. 



^ R 



r/2. In this limit Eq. (43) simplifies to 



St = 



JF^ 



n±[(^ 



/l — Jr 
-e±6)2H 



(F/2)^ 



(58) 



More generally, the sign of st is determined by the con- 
dition 



sgn[st] = sgn [(FtF^ - Flrj,) (/^ - fn) 



(59) 



Thus, when we consider a f-polarized left lead with 
rl^fl > rj^r]^, a current fiowing from left to right, 
fJ-L > fJ-R, results in S( < and thus antialignment of 
magnetic moment and magnetic field. For the opposite 
spin-polarization, the spin-torque tends to align the mag- 
netic moment with the magnetic field, as sketched in Fig. 

m 

For a given spin-polarization, inverting the direction 
of the current can switch the orientation of the magnetic 
moment in the same way. This is studied by solving nu- 
merically the equation of motion for the molecular spin 
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in the strictly adiabatic limit, hence neglecting Gilbert 
damping and fluctuations, in the presence of strongly 
polarized leads. In Fig. |8] we show the time evolution 
of the molecular spin initially slightly deviating from the 
magnetic field axis for two different voltages. Clearly, the 
motion of the molecular spin is determined by the direc- 
tion of the current through the molecule, showing that 
inverting the bias voltage causes spin-flips in this setup. 



VIII. SUMMARY AND CONCLUSIONS 

In this work we have considered an anisotropic mag- 
netic molecule in a single-molecule junction in which 
conduction electrons couple via exchange to the local- 
ized magnetic moment. The resulting current-induced 
torques have been analyzed by means of the non- 
equilibrium Born-Oppenheimer approximation, which 
gives rise to Langevin dynamics of the magnetic moment, 
described by a generalized Landau-Lifshitz-Gilbert equa- 
tion. This approximation is valid in the high-current 
limit when the precessional frequency of the molecular 
spin is small compared to the electronic time scales. Un- 
like previous works, our approach does not follow a per- 
turbative route either in the tunneling between leads and 
the molecule or in the coupling between the electronic 
spin and the molecular magnetic moment. Accordingly, 
we can render the full dependence of the parameters of 
the LLG equation on the state of the molecular moment 
as well as on the applied bias and gate voltages. 

The strictly adiabatic approximation causes a mean 
torque exerted by the conduction electrons, while retar- 
dation effects result in a renormalization of the precession 
frequency and Gilbert damping. In addition, equilibrium 
and non-equilibrium fluctuations of the current cause a 
fluctuating (Langevin) torque. We have expressed these 
torques in terms of the electronic Green's functions and 
have related them to scattering theory, in the latter case 
extending earlier work to include an applied bias voltage. 
We have concluded that in general out-of-equilibrium sit- 
uations the conduction electrons can transfer energy to 
the localized moment by the fluctuations and, in the pres- 
ence of spin-polarized leads, via a non-conservative (spin- 
transfer) torque and/or negative damping. 

These mechanisms allow one to use the anisotropic 
magnetic molecule in an external magnetic field as a 
molecular switch which can be read out via the back- 
action of the molecular spin on the transport current. 
When the molecule is attached to metallic leads in a 
uniaxial setup, we have turned the Langevin equation 
into a Fokker-Planck equation allowing us to calculate 
the switching rates between the two stable spin orienta- 
tions. Transitions between these states are driven by the 
fluctuations which we have analyzed-^ in addition to the 
mean torque, damping, and the current- as a function 
of the applied gate and bias voltages and the orientation 
of the molecular spin. In the presence of spin-polarized 
leads, the switching dynamics is dominated by the non- 



conservative (spin-transfer) part of the current-induced 
torque, which enables switching between the spin orien- 
tations by reversing the direction of the electronic cur- 
rent. 

The above mentioned features of the dynamics of the 
local magnetic moment are also common in layered spin- 
tronic devices. However, in the present case, the differ- 
ent coefficients that govern the dynamics of the molecu- 
lar magnetic moment show a strong dependence on the 
bias voltage determined by the electronic structure of the 
molecule (see Fig. 3). The latter property also determines 
the behavior of the electronic current, where features of 
the dynamics of the magnetic moment take place in com- 
bination with coherent tunneling of molecular systems, as 
signalized for instance in the current and the differential 
conductance (see Fig. [6| . 

We have considered a generic and standard model for 
the molecule which applies to a wide type of molecular 
systems, provided that a sufficiently large current flows 
through the molecule and that the magnetic moment is 
sufficiently large to fulfill the adiabatic condition assumed 
in the NEBO treatment. In particular, good candidates 
can be the Mni2- or Fcg-based devices. These systems 
are described by microscopic Hamiltonians of the type we 
considered in this work, and have rigid magnetic cores 
with M = 10, and an anisotropy barrier DM^ of the 
order of a few meV.^ Classical descriptions of their 
spin dynamics have been presented for these molecules 
in contact to phononic environments.'^ The crucial pa- 
rameters in order to achieve the adiabatic regime in our 
setup, should be a good enough contact to the electrodes 
and a sufficiently high applied bias voltage, leading to 
a short dwell time of the electrons in the molecule. To 
be more specific, we estimate for the Mni2- or Feg sys- 
tems with a rather large magnetic anisotropy that the 
Born-Oppenheimer approximation can be applied when 
the current through the device exceeds ^ lOnA. 
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Appendix A: Green's functions, s°(i,M), 'y'J'{t,M) 
and noise correlator 

We approximate the self-energy to be independent of 
energy. In this wide band limit Eq. ( 13 1 becomes 



I]f,,(a;)^Re(sf_,) 



7^^ ly, 



(Al) 
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with the approximately constant density of states resulting in Eq. ( 40 1 in the case of axial symmetry. The 



Va,a(!-o) ^ and \w\ = 
troduce the abbreviations 



(ki 



WflH /2. We in- 



= 6o+^Re(l]«J (A2) 

and we will use the notation ~ [Ej^ ± S|*]/2 and 

^c,s = Ea ^a,c(s) , with r^^^cCs) = 0^0,^ ± taking 
into account possibly spin-polarized leads. 

From Eq. (18) we find for the frozen retarded Green's 
function 



1 



G«(M,u;) = -^ 



w 1 + b • cr 
ll + cr-b/6 ll-cr-b/& 



-b 



uj + b 



(A4) 



with u) = (lu — €o — S^) ~ w— e+irc/2. Here, we include 
the antisymmetric part of the self-energy in the effective 
magnetic field, 



b(t)-i(JM(t)+5,B) 



(A5) 



After some algebra we find the following expression for 
the lesser Green's function (|2Tl: 



(M, w) =Gf (M, w) + G< (M, w)b • cr + G< (M, w)cr" 



+ G<(M,cj)(T- (e^ X b), 
where the coefficients are given by 
1 



(A6) 



G<(M,^) 
G<{M,cj) 

G<{M,uj)=- 



|A(M,c.)|2 



~,|2 



S]<(w)2[Re[w]6;, 



-]}, 



:^^^^{E<HRe[^]+I]<M54, 

1 r r 



|A(M,c.)|2 
+ S<H[|<ip-|bp 



2 



A(M,a.)| 



js<Hr,-i]<Hrj. (A?) 



We use A(M, uj) ^ - iP a.nd h c::^ h ~ i{Ts/2)e^. Sub- 
stituting the above expressions for T^s, it can be seen 
that Gf (M,a;) = 0, for r^^o- = ^R,a, implying that this 
component of the Green's function contributes only for 
polarized leads. Note that the corresponding expressions 
for the larger Green's function G5'(M,a;) are obtained 
by replacing S<^ by E>^ = (1 - fa{i^))^a,c(s) in 

the expressions above. 

Using the Green's functions expressions, we find for 
the mean value of the electronic spin at the molecule 



s"(M) 



dw 
2Td 



{G<b + G<e, + G<(e, x b)}, (A8) 



explicit expression for the component parallel to M reads 
iJ f duj 



sm{M) 



M / 27r 



G<{M,oj). 



(A9) 



The correction due to retardation effects associated 
with the slow dynamics of the molecular spin are cap- 
tured by the matrix 7'*^, see Eq. (24). The symmetric 
part of this matrix. 



7f(M) 



J2 

4 



duj 
2n 



tr[a'a^G>a'=G<i 



(AlO) 



describes Gilbert damping of the molecular spin, induced 
by the coupling to the electrons. The antisymmetric part 
of the matrix 7 is given by 



J2 

4 



duj 

271 



tiyd^ (G« + G'o 



<y''G<] 



(All) 



Considering a setup with unpolarized leads and the ex- 
ternal magnetic field pointing along the anisotropy axis, 
hence b = ( JM + geBe^) /2, Eq. (|A10[) becomes 



7f(M) 



" 4 
' 2 



r„r^(-9.//3) 



duj X ^ 

^t^U± [(^-6±6)2 + (r/2)^ 

duj 

a/3 



[(l^ - e ± 6)2 + (r/2)2 



>bkbi. 



(A12) 



This will be decomposed into a term proportional to the 
unit matrix and a projector onto the z-axis, as described 
in Sec. |v] Note that the sign of the eigenvalues of 7^' 
is fixed, corresponding to damping in and out-of equi- 
librium. As described in the main text, the pref actor a 
in Eq. (ISll is given by a{M) ^ 1/{1 + C- M), with C 



defined in Eq. ( 29 1 . This becomes 

J2r2r«/a(c^-e) 



ciM)= fpy- 



[(c^ - e ± 6)2 + (r/2): 



jib, 
(A13) 



where we have inserted Gq , Eq. (A6I, and the corre- 



sponding expression for Gq into Eq. (All). 



Appendix B: Fokker-Planck equation 

In this appendix we derive the Fokker-Planck equa- 
tion from the Langevin equation and obtain an expres- 
sion for the mean first passage time, following standard 
arguments.!^ 

We note that the probability distribution for the molec- 
ular spin is conserved for all t, J dM /(M, t) — 1. Hence, 
we can write a continuity equation for the probability dis- 
tribution, 

5t/(M,t) + aM- (M/(M,t)) =0. (Bl) 
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Inserting Eq. (31 1 for M we get 



dtf{M, t) = -L /(M, t)~dM- (i)/(M, t)) , (B2) 

where ^{t) = M x SB and the differential operator L is 
defined via its action on the function /(M, t) as 



/ 



(B3) 



Lf ^Om- (aM X 
From this follows the implicit solution 
f{M,t) =e-*^/(M,0) 

dt'e-(*-*')^9M-(l(t')/(M,t')). (B4) 



Inserting this again in Eq. (B2) and averaging over noise, 
denoted by P{M,t) = (j{M,t)), yields the Fokker- 
Planck equation 

dtP{M, t) = -L F(M, t) + • (a'D) • d^P{m, t) 



FP{M,t), 



(B5) 



where we use that the noise is Gaussian and delta- 
function correlated, (Cfc(06(*')) = D^^5{t — t') and in- 
troduce the Fokker-Planck operator F. 



We consider the distribution P(M, t) of M which have 
been at Mq at time t = Q and are inside a given volume 
at time t. The mean first passage time is then given by 



r(Mo) = j dtt j dm 



-dP{m,t) 
Jt 



(B6) 



where — / dM ^'^^J^'*'' is the distribution of first passage 
times and / dMP(M, t) gives the number of M which 
are still in the volume of consideration at time t. The 
distribution of M is P[M,t) = e^^ 5{M. - Mq) with 
P(M,t) = when M is at t he boundary of the volume. 
We insert this into Eq. ( |B6[ ) so that after integration by 
parts 

r(Mo)= j dt j Me*^(5(M - Mo) = j dte'^h, (B7) 

with the adjoint Fokker Planck operator . This results 
in the differential equation 



i^V(M) = -1 



(B8) 



for the mean first passage time with an absorbing bound- 
ary condition. 
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